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Abstract 

The Noether symmetry of a generic f{R) cosmological model is investigated by utilizing the 
behavior of the corresponding Lagrangian under the infinitesimal generators of the desired sym- 
metry. We expHcitly calculate the form of f{R) for which such symmetries exist. It is shown that 
the resulting form of f{R) yields a power law expansion for the cosmological scale factor. We also 
obtain the effective equation of state parameter for the corresponding cosmology and show that 
our model can provide a gravitational alternative to the quintessence. 
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^ ; 1 Introduction 

>, 

On , Observations of type la supernovae [1] and cosmic microwave background [2] have revealed that the 

universe is undergoing an accelerating phase in its expansion. The explanation of this acceleration 
fT^ ' in the context of general relativity has stimulated a myriad of ideas, the most notable of which is 

.^ . the introduction of a mysterious cosmic fluid, the so-called dark energy [3]. In more recent times, 

C^ I modified theories of gravity, constructed by adding correction terms to the usual Einstein-Hilbert 

action, have opened a new window to study the accelerated expansion of the universe where it has 
been shown that such correction terms could give rise to accelerating solutions of the field equations 
without having to invoke concepts such as dark energy [4]. 
p\ ' In a more general setting, one can use a generic function, f{R), instead of the usual Ricci scalar R 

S . as the action for the model. Such f{R) gravity theories have been extensively studied in the literature 

over the past few years. In finding the dynamical equations of motion one can vary the action 
with respect to the metric (metric formalism), or view the metric and connections as independent 
dynamical variables and vary the action with respect to both independently (Palatini formalism) 
[5]. In this theory, the Palatini form of the action is shown to be equivalent to a scalar-tensor type 
theory from which the scalar field kinetic energy is absent. This is achieved by introducing a conformal 
transformation in which the conformal factor is taken as an auxiliary scalar field [6]. As is well known, 
in the usual Einstein-Hilbert action these two approaches give the same field equations. However, in 
f{R) gravity the Palatini formalism leads to different dynamical equations due to nonlinear terms 
in the action. There is also a third version of f{R) gravity in which the Lagrangian of the matter 
depends on the connections of the metric (metric-affine formalism) [7]. 

In this letter we consider an (n-l-l)-dimensional flat FRW space-time in the framework of the 
metric formalism of f{R) gravity. Following [8], we set up an effective Lagrangian in which the scale 
factor a and Ricci scalar R play the role of independent dynamical variables. This Lagrangian is 
constructed in such a way that its variation with respect to a and R yields the correct equations of 
motion as that of an action with a generic f{R) mentioned above. The form of the function f{R) 
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appearing in the modified action is then found by demanding that the Lagrangian admits the desired 
Noether symmetry [9]. By the Noether symmetry of a given minisuperspace cosmological model we 
mean that there exists a vector field X, as the infinitesimal generator of the symmetry on the tangent 
space of the configuration space such that the Lie derivative of the Lagrangian with respect to this 
vector field vanishes. For a study of the Noether symmetry in various cosmological models see [10]. 
We shall see that by demanding the Noether symmetry as a feature of the Lagrangian of the model 
under consideration, we can obtain the explicit form of the function f{R)- Since the existence of a 
symmetry results in a constant of motion, we can integrate the field equations which would then lead 
to a power law expansion of the universe. 

2 The phase space of the model 

In this section we consider a spatially flat FRW cosmology within the framework of f{R) gravity. 
Since our goal is to study models which exhibit Noether symmetry, we do not include any matter 
contribution in the action. Let us start from the (n+l)-dimensional action (we work in units where 
c = IGttG = 1) 

S= I <P+^x./^gf{R), (1) 



where R is the scalar curvature and f{R) is an arbitrary function of R. By varying the above action 
with respect to metric we obtain the equation of motion as 

\g^.f{R) - R^uf'iR) + V^V,/'(i?) - g^uOf'iR) = 0, (2) 

where a prime represents differentiation with respect to R. We assume that the geometry of space- 
time is described by the flat FRW metric which seems to be consistent with the present cosmological 
observations 

n 

ds'^ = -df + a\t)Y,idx'f- (3) 

With this background geometry the fleld equations read 



in-l)H^^^H^ = -^ 



f"'R2 + („ _ i^HRf" + f'R + i(/ - Rf) 



(4) 



H' 



— ^— - [ifR - /) - 2nRHf"] , (5) 

n(n — l)j' L J 

where H = a/a is the Hubble parameter and a dot represents differentiation with respect to t. To 
study the symmetries of the minisuperspace under consideration, we need an effective Lagrangian 
for the model whose variation with respect to its dynamical variables yields the correct equations of 
motion. Following [8], by considering the action described above as representing a dynamical system 
in which the scale factor a and scalar curvature R play the role of independent dynamical variables, 
we can rewrite action (1) as [11] 



S= f dtC{a, h, R,R)= I dt\ a"/(i?) - A 



• 2 

a a 

R — n(n — 1)—^ — 2n— 

a^ a 



(6) 



where we introduce the deflnition of R in terms of a and its derivatives as a constraint. This procedure 
allows us to remove the second order derivatives from action (6). The Lagrange multiplier A can be 
obtained by variation with respect to R, that is, A = a"'f'{R). Thus, we obtain the following 
Lagrangian for the model [8, 11] 

C{a, a, R, R) = n(n - l)a^a''~^f' + 2naRa''-^ f" + a'\f R - /). (7) 



The momenta conjugate to variables a and R are 

Pa = ^ = 2n(n - l)da"-2/' + 2na''-^Rf", (8) 

BC 

PR = -5 = 2nda"-V". (9) 

oR 

The Hamiltonian corresponding to Lagrangian (7) can then be written in terms of a, a, R and R as 

nia, a, R, R) = n{n - l)a'^a'^-^f' + 2ndi?a"- V" - a"(/'i? - /). (10) 

Therefore, our cosmological setting is equivalent to a dynamical system where the phase space is 
spanned by {a,R,pa,PR} with Lagrangian (7) describing the dynamics with respect to time t. Now, 
it is easy to see that variation of Lagrangian (7) with respect to R gives the following well-known 
relation for the scalar curvature 

• 2 

R = n(n-l)% + 2n-, (11) 

a^ a 

while variation with respect to a yields the field equation (4). Also, equation (5) is nothing but the 
zero energy condition Ti. = (Hamiltonian constraint). A quantum f{R) cosmological model based 
on Hamiltonian (10) can be found in [12]. 

The setup for constructing the phase space and writing the Lagrangian is now complete. In the 
next section we shall use the above Lagrangian to find the form of f{R) that would admit a Noether 
symmetry. 

3 The Noether symmetry 

As is well known, Noether symmetry approach is a powerful tool in finding the solution to a given 
Lagrangian, including the one presented above. In this approach, one is concerned with finding 
the cyclic variables related to conserved quantities and consequently reducing the dynamics of the 
system to a manageable one. The investigation of Noether symmetry in the model presented above is 
therefore the goal we shall pursue in this section. Here our aim is to find the function f{R) such that 
the corresponding Lagrangian exhibits the desired symmetry. Following [9], we define the Noether 

symmetry induced on the model by a vector field X on the tangent space TQ = (a,R,a,R] of the 
configuration space Q = (a, R) of Lagrangian (7) through 

^^ d ^ d da d dd d , , 

such that the Lie derivative of the Lagrangian with respect to this vector field vanishes 

LxC = 0. (13) 

In (12), a and /3 are functions of a and R and ^ represents the Lie derivative along the dynamical 

vector field, that is, 

d . d ■ d ,, ^, 

Jt='0-a^''0R- ^''^ 

It is easy to find the constants of motion corresponding to such a symmetry. Indeed, equation (13) 

can be rewritten as 

^ ^ f dc dadc\ f ac d/3dc\ ^ 



Noting that ^ = -^, we have 



dpa da \ ( r,dPR dp 



which yields 



Thus, the constants of motion are 



dt 



{apa + Ppr) = 0. 



Q = apa + Ppr. 



(17) 



(18) 



In order to obtain the functions a and /? we use equation (15). In general this equation gives a 
quadratic polynomial in terms of d and R with coefficients being partial derivatives of a and (3 with 
respect to the configuration variables a and R. Thus, the resulting expression is identically equal to 
zero if and only if these coefficients are zero. This leads to a system of partial differential equations 
for Q and /3. For Lagrangian (7), condition (15) results in 



nd'-^ifR - f)a + a"/"i?/? + 2nR^ ( a"" V — ^ 



+nd 

+2nRa 
= 0, 



dR) 



,n-l f//f, 



n-1 ftl 



da 



da 

■a-2 r; da 



da 



(n - l)a"-^/"a + a^-^'"/? + «"" V'V + (^ " !)«""/ 7^5 + ^ flHS 



da 



dR 



dR 



(19) 



which leads to the following system of equations 



n{f'R - f)a + af'Rp 

f,fda 
■' dR 

(n - l)(n - 2) fa + (n - l)af"(3 + 2(n - l)a/'|^ + 2aV"^ 
(n-l)ra + ar/. + ar|^ + (n-l)/|| + ar|| 

Equation (21) gives /" = or 4^ = 0. In the case where /" = 0, we obtain 

fiR) = ciR + C2. 



0, 
0, 

0, 

0. 



(20) 

(21) 

(22) 
(23) 

(24) 



Substituting this result into equation (20) we get C2 = and recover the usual Einstein-Hilbert gravity 
without cosmological constant for which f{R) = R. Also, since equation (23) results in 4^ = 0, from 
equation (22) we obtain the following differential equation for a{a) 



{n-2)a{a) + 2a-^ = 0, 
da 



with solution 



a(a} 



n — 2 

a 2 . 



(25) 



(26) 



Therefore, for any arbitrary function /3(a, R), the following vector field represents a Noether symmetry 
for the minisuperspace of the flat FRW cosmology in Einstein-Hilbert gravity 



n-2 d n — 2._n_d ^ d ■ d 



(27) 



Since from (9) we have pR = 0, the corresponding constant of motion is Q = apa, which yields 

a{t) ~ t^/n, (28) 



a '^ aa = const. 



Although, this solution for the scale factor satisfies the field equation (4), it does not satisfy the 
Hamiltonian constraint (5). This is not surprising since it is well known that the flat FRW cosmology 



with Einstein-Hilbert action (without the cosmological constant) has no vacuum solution. Thus, we 
remove the case /" = from our consideration. When /" 7^ 0, from equation (21) we obtain 



da 

an-"- P^' 



Under this condition equation (20) results in 



«(«) = ^f^^/3K«)- (30) 



If one uses this expression in equation (22) to eliminate a{a), one obtains 

d/3 n{n - 1)/ 



da 2a{f'R - nf) 



P{a,R). (31) 



To solve this equation we assume that the function (3{a, R) can be written in the form f3{a, R) = 
A{a)B{R), where A and B are functions only of a and R respectively. Substituting this ansatz for 
P{a,R) into equation (31), we obtain 

2adA ^ n(n-l)/ 

Ada f'R-nf ^ ' 

Since the left-hand side of this equation is a function of a only while the right-hand side is a function 

oi R, we should have 

n{n- 1)/ 



f'R - nf 
which results in 



c = Const., (33) 



n(n-\-c—l) 

fiR) = R^^, (34) 

where the constant c 7^ —n is required to have /" 7^ 0. On the other hand equation (32), with its 
right hand-side equal to c, has the solution 

'^adA ^. . „/2 ,„^, 

— = c^ A(a) = a"'. (35) 

a da 

Now, using (34) and P{a,R) = a'^'^B{R) in equation (30), we find 

a{a) = _ " + c-l ^-ia§+i^(^). (36) 

c 

Since a{a) should be a function of a only, from the above expression for a{a) we conclude that 
B{R) = R and thus obtain 

p{a,R) = Ra^/'^, a{a) = _ ^ + g-l ^f+i_ (37) 

c 

To determine the constant c, we note that relations (34) and (37) should satisfy equation (23). Thus, 
by substituting these results into equation (23) we find 

c=-l-n, (38) 

which completes our solutions as 

Zi n—1 , , n + 1 

a{a) = -a —, p{a,R) = Ra —, (39) 

and 

2n 

f{R)=R^+T. (40) 



Therefore, in the context of f{R) cosmology, a flat FRW metric has Noether symmetry if the corre- 
sponding action is given by equation (40) and the Noether symmetry is generated by the following 
vector field 



2 _ii^ d „ _2i±i d n — 1 _ii±i d 
n + 1 oa on n + 1 oa 

„n+l n + 1^. _n+3\ d 

+ [ Ra 2 Raa 2 — ^. 

2 J dR 



(41) 



To obtain the corresponding cosmology resulting from this type of /(R), we note that the existence 
of Noether symmetry (41) implies the existence of a constant of motion Q = apa + Ppr- Hence, using 
equation (8) and (9) we have 



Q 



4n^(n — 1). n^ li^ 8n^(n — 1) n-i 



(n + l)2 
which may be rewritten in the form 

Q = 



-da 2 i?"+i 



{n + lf 



a 2 RR n+i 



8n? d f n-i "-i\ 

a 2 i?n+l ^ 



(n + 1)2 dt 

and can be immediately integrated with the result, assuming a{t = 0) = 

(n+l)2 



n— 1 n—1 

a 2 i?n+l 



8n2 



-Qt. 



On the other hand, the Hamiltonian constraint 7i = with 7i given by (10) gives 

4n2 



2n^a'^R-^ + 



n+ 1 



-aRaR = a 



To obtain the scale factor from the above relation, note that equation (44) yields 



R 



(n+1) 
8n2 



-Q 



n + l 

"""^ ^n + 1 _2_ _n+l n+l H+l n+3 

t"-ia 2 t"-ida 2 

n-1 2 



which, upon substitution into relation (45), gives 



n — 1 
where, after integration we obtain 



(n + l) 
8n2 



-Q 



n + l 
' n-1 



n-1 

da 2 



2n 



a(t) ~ t 



2(3n.-l) 
n2-l 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



(48) 



Therefore, in the context of our f{R) = R"+^ gravity, the universe evolves with a power law expansion 
(note that for any n > 1 the power of the scale factor is positive). It is remarkable from (48) that 
the condition under which the universe would accelerate is ■? , > 1, that is, n < 5. This means 
that models with spatial dimension n < 5 obey an accelerated power law expansion while for n > 5 
a decelerated expansion occurs. 



4 The equation of state parameter 

One of the advantages of f{R) theories of gravity is that they can describe and provide gravitational 
alternative for dark energy. Indeed, the equations of motion resulting from f{R) gravity admit such 



solutions which predict the same accelerated expansion as those resulting from the usual Einstein- 
Hilbert gravity with dark energy. In other words, in the context of f{R) cosmology one can describe 
the accelerated expansion of the universe without the necessity to introduce the exotic fluid with a 
negative equation of state (EoS) parameter w. To introduce an effective EoS parameter in our model, 
one may compare the field equations (4) and (5) with the usual field equations of a flat FRW universe 
filled with a perfect fluid with EoS P = wp, that is, with 3H^ ~ p and -2H - 3H'^ ~ P. This 
comparison shows that in our model one may introduce the effective energy density and pressure as 
follows 

P=\{Rf'-f)-nHRf'\ (49) 

P = f'R^ + (n - l)HRf" + Rf" + \{f- Rf) . (50) 

Therefore, we can define the effective EoS parameter w^jf as 

P f"'R' + {n- l)HRf" + Rf" + \U- Rf) 

Weff = — = q ■■ . (51) 

" P \{Rf'-f)-nHRf" 

To transform the above expression to a more manageable form, consider the quantity J- = f'{R) [8]. 
In terms of T, equation (51) takes the form 

2 T-HT 
'■' n{n — 1) J^H^ 

For the problem at hand, taking into account that 

f{R) = ^R^, (53) 

n + 1 

and with the help of equation (44) we obtain 

72+1 on — 1 

T=-^^Qt-^—. (54) 

4n 

Thus, substituting (54) and H = aja from (48) we are led to the following EoS parameter 

(n-l)(n + l) , , 

n(3n — Ij 

Now, it is easy to see that — 1 < w^fj < 0, which is the characteristic of one type of dark energy, 
the so-called quintessence. Note that when w^ff is less than —1, the EoS describes another type 
of dark energy known as phantom. As is clear from equation (55), the effective EoS parameter is 
always greater than —1 and thus in the model under consideration here, the phantom phase cannot 
be accounted for. 

5 Conclusions 

In this letter we have studied a generic f{R) cosmological model with an eye to Noether symmetry. 
For the background geometry, we have considered a flat (n+l)-dimensional FRW metric and derived 
the general equations of motion in this background. The phase space was then constructed by taking 
the scale factor a and Ricci scalar R as the independent dynamical variables. The Lagrangian of 
the model in the configuration space spanned by {a, R} is so constructed such that its variation 
with respect to these dynamical variables yields the correct field equations. The existence of Noether 
symmetry implies that the Lie derivative of this Lagrangian with respect to the infinitesimal generator 
of the desired symmetry vanishes. By applying this condition to the Lagrangian of the model, we have 



obtained the explicit form of the corresponding f{R) function. We have shown that this form of f{R) 
results in a power law expansion for the scale factor of the universe and the expansion is accelerating 
in the case of n < 5. We have also presented an effective EoS parameter for our f{R) cosmology 
model. Our analysis shows that the EoS parameter is always greater than —1 and thus the so-called 
phantom dark energy cannot be described in this kind of modified gravity. On the other hand, the 
EoS parameter is restricted to the interval — 1 < Wgff < 0. Since this interval is characteristic of the 
quintessence dark energy, our f{R) model can provide a natural gravitational alternative for this kind 
of dark energy without the necessity to introduce an exotic fluid with a negative EoS parameter. 

Acknowledgement 

The author would like to thank H.R. Sepangi and N. Khosravi for a careful reading of the manuscript 
and useful comments. 

References 

[1] A. G. Riess et. al, Astron. J. 116 (1998) 1009 (arXiv: astro-ph/9805201) 
S. Perlmutter et. al, Astrophys. J. 517 (1999) 565 (arXiv: astro-ph/9812133) 
D. N. Spergel et. al, Astrophys. J. Suppl. 148 (2003) 175 (arXiv: astro-ph/0302209) 

[2] C. L. Bennett et. al, Astrophys. J. Suppl. 148 (2003) 1 (arXiv: astro-ph/0302207) 
C. B. Netterfield et. al, Astrophys. J. 571 (2002) 604 (arXiv: astro-ph/0104460) 
N. W. Halverson et. al, Astrophys. J. 568 (2002) 38 (arXiv: astro-ph/0104489) 

[3] S. M. Carroh, Living Rev. Rel. 4 (2001) 1 (astro-ph/0004075) 

T. Padmanabhan, Phys. Rep. 380 (2003) 235 (arXiv: hep-th/0212290) 

M. Tegmark et. al, Phys. Rev. D 69 (2004) 103501 (arXiv: astro-ph/0310723) 

V. Faraoni, Class. Quantum Grav. 22 (2005) 3235 (arXiv: gr-qc/0506095) 

S. Capozziello, S. Nojiri and S. D. Odintsov, Phys. Lett. B 632 (2006) 597 (arXiv: hep- 

th/0507182) 

S. Nojiri and S. D. Odintsov, Gen. Rel. Grav. 38 (2006) 1285 (arXiv: hep-th/0506212) 

S. Nojiri, S. D. Odintsov and S. Tsujikawa, Phys. Rev. D 71 (2005) 063004 (arXiv: hep- 

th/0501025) 

M. Amarzguioui, O. Elgaroy, D. F. Mota and T. Multamaki, Astron. Astrophys. 454 (2006) 707 

(arXiv: astro-ph/0510519) 

S. Nojiri and S. D. Odintsov, Int. J. Geom. Meth. Mod. Phys. 4 (2007) 115 (arXiv: hep- 

th/0601213) 

[4] C. Deffayet, Phys. Lett. B 502 (2001) 199 (arXiv: hep-th/0010186) 

J. S. Alcaniz, Phys. Rev. D 65 (2002) 123514 (arXiv: astro-ph/0202492) 

S. M. Carroh, V. Duvvuri, M. Trodden and M. Turner, Phys. Rev. D 70 (2004) 043528 (arXiv: 
astro-ph/0306438) 

S. Nojiri and S. D. Odintsov, Phys. Lett. B 576 (2003) 5 (arXiv: hep-th/0307071) 
K. Atazadeh and H. R. Sepangi, Int. J. Mod. Phys. D 16 (2007) 687 (arXiv: gr-qc/0602028) 
S. Nojiri, S. D. Odintsov and M. Sami, Phys. Rev. D 74 (2006) 046004 (arXiv: hep-th/0605039) 
I. Navarro and K. Van Acoleyen, J. Cosmol. Astropart. Phys. JCAP 0603 (2006) 008 (arXiv: 
gr-qc/0511045) 

S. Capozziello, Int. J. Mod. Phys. D 11 (2002) 483 (arXiv: gr-qc/0201033) 
S. Capozziello, S. Carloni and A. Troisi, Recent Research Developments in Astronomy and As- 
trophysics RSP/AA/21-2003 (arXiv: astro-ph/0303041) 

S. Nojiri and S. D. Odintsov, Phys. Rev. D 68 (2003) 123512 (arXiv:hep-th/0307288) 
S. Nojiri and S. D. Odintsov, (arXiv: 0801.4843 [astro-ph]) 



[5] D. N. Vollick, Phys. Rev. D 68 (2003) 063510 (arXiv: astro-ph/0306630) 

X. Meng and P. Wang, Class. Quantum Grav. 20 (2003) 4949 (arXiv: astro-ph/0307354) 

X. Meng and P. Wang, Palatini formulation of modified gravity with squared scalar curvature 

(arXiv: astro-ph/0308284) 

[6] E. E. Flanagan, Class. Quantum Crav. 21 (2003) 417 (arXiv: gr-qc/0309015) 
E. E. Flanagan, Phys. Rev. Lett. 92 (2004) 071101 (arXix: astro-ph/0308111) 
G. Magnano, L. M. Sokolowski, Phys. Rev. D 50 (1994) 5039 (arXiv: gr-qc/9312008) 

[7] T. P. Sotiriou and S. Liberati, Ann. Phys. 322 (2007) 935 (arXiv: gr-qc/0604006) 
N. J. Poplawski, Class. Quantum Crav. 23 (2006) 2011 (arXiv: gr-qc/0510007) 

[8] J. C. C. de Souza and V. Faraoni, Class. Quantum Crav. 24 (2007) 3637 (arXiv: 0706.1223 
[gr-qc]) 

[9] M. Demianski, R. de Ritis, C. Rubano and P. Scudellaro, Phys. Rev. D 46 (1992) 1391 

S. Capozziello, G. Marmo, C. Rubano and P. Scudellaro, Int. .J. Mod. Phys. D 6 (1997) 491 
(arXiv: gr-qc/9606050) 

S. Capozziello, S. Nesseris and L. Perivolaropoulos, J. Cosmol. Astropart. Phys. JCAP 0712 
(2007) 009 (arXiv:0705.3586 [astro-ph]) 

S. Capozziello and A. De Felice, f{R) cosmology by Noether's symmetry (arXiv: 0804.2163 [gr- 
qc]) 

[10] S. Capozziello, V. I. Man'ko, G. Marmo and C. Stornaiolo, Tomographic Representation of 
Minisuperspace Quantum Cosmology and Noether Symmetries (arXiv: 0706.3018 [gr-qc]) 
S. Capozziello, A. Stabile and A. Troisi, Class. Quantum Crav. 24 (2007) 2153 (arXiv: gr- 
qc/0703067) 

A. K. Sanyal, B. Modak, C. Rubano and E. Piedipalumbo, Cen. Rel. Grav. 37 (2005) 407 (arXiv: 
astro-ph/0310610) 

A. K. Sanyal, Phys. Lett. B 524 (2002) 177 (arXiv: gr-qc/0107053) 
A. K. Sanyal and B. M. Modak, Class. Quantum Grav. 18 (2001) 3767 (arXiv: gr-qc/0107052) 

A. Bonanno, G. Esposito, C. Rubano and P. Scudellaro, Gen. Rel. Crav. 39 (2007) 189 (arXiv: 
astro-ph/0612091) 

B. Vakili, N. Khosravi and H. R. Sepangi, Class. Quantum Crav. 24 (2007) 931 (arXiv: gr- 
qc/0701075) 

[11] S. Capozziello and G. Lambiase, Gen. Rel. Crav. 32 (2000) 673 (arXiv: gr-qc/9912083) 

[12] A. Shojai and F. Shojai, f(R) Quantum Cosmology, (arXiv: 0801.3496 [gr-qc]) 



